GLUONS IN THE LATTICE SU(2) CLASSICAL GAUGE FIELD 



YING CHEN, JI-MIN WU, BING HE 

Institute of High Energy Physics, Chinese Academy of Sciences, 
P.O.Box 918-4, Beijing 100039, P. R. China 
E-mail: cheny@hptc5.ihep.ac.cn 

The SU(2) gluonic correlation functions, glueball effective masses in the J p = + , 
2+ and 0~ channels were calculated from the lattice classical gauge configurations 
which were obtained by smoothing the thermal gauge configurations through the 
improved cooling method. The instanton-induced attractive force in the + chan- 
nel and the repulsive force in the — channel are confirmed in the Monte Carlo 
simulation. There is evidence that the instanton vacuum contribution to the + 
glueball mass is significant. 



1 Introduction 

QCD predicts there exist glueballs. The structure and formation mechanism 
of glueballs are so extricated that it must be explored nonperturbatively. Ex- 
tensive Monte Carlo simulations have been performed and a few statements 
seem to be firmly established in the glueball sectonl: i) the scalar glueball is the 
lightest with the mass in the 1.6-1.8 GeV range; ii) the tensor glueball and the 
pseudoscalar one are significantly heavier than the scalar one with, the mass 
ratio 77i 2 ++ /m Q ++ ~ 1.4 and to -+/to ++ ~ 1.5— 1.8, respectiveljfl; iii) The 
scalar is much smaller than other glueballs (r ++ ~ 0.2/m, r 2 ++ ~ 0.8/rrH). 
From these statements one can find that glueballs are much heavier than typ- 
ical quark-model hadrons and the size difference betweensi the scalar and the 
tensor glueball is much larger than that of the similar mesons( a similar mea- 
surement for the 7r and p mesons gives 0.32/m and 0.45/m). This indicates 
that the spin-dependent forces between gluons are stronger than that between 
quarks. 

The information about the gluonic interactions can be obtained by study- 
ing the correlation functions oLgluonic operators with the relevant quantum 
numbers. Schafer and ShuryalH calculated the correlation functions and the 
Bethe Salpeter amplitudes in an instanton-based model of the QCD vacuum. 
Their results show that instantons lead to a strong attractive force in the 
J PC = ++ channel, thus the scalar glueball is much smaller than other glue- 
balls. In the h channel the corresponding force is repulsive, and in the 2 ++ 
case it is absent. 

Recently the cooling method, a smoothing algorithm for the gauge config- 
urations, has been developed to extract classical contents( such as instantons) 
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of QCD vacuum on the lattice. This motivates us to study the relationship 
of instantons and glueballs by the lattice simulation. The goal of this work is 
to study the instanton effects in the glueball masses and gluonic correlation 
functions in the lattice SU(2) gauge theory. The tadpole improved Symanzik's 
action] is used on an anisotropic lattice to perform the Monte Carlo simu- 
lation ^s accurate as possible on coarse lattices. The smearing and fuzzy 
schemtH are applied to the glueball operators to increase the signal-to-noise 
ratio. These schemes^have been verified to be successful in Morningstar and 
Peardon's calculations and in our work. We employ the improved cooling to 
extract the instanton configurations. The gluonic correlators are measured 
during the cooling procedure and the glueball masses are extracted by the ex- 
ponential fall-off of these correlators with respect to the Euclidean time. To 
investigate the spin dependence of the instanton effects, the behaviors of the 
gluonic vacuum correlation functions are also explored in this work. 

The simulation details, the results and the discussion are given in section 
II. Section III is the conclusion and summary. 



2 Lattice Calculation 

If the lattice action is properly improved, the costs of the Monte Carlo simula- 
tion will be greatly reduced and the numerical study can be performed on much 
coarser lattice system. In this work, the gluonic correlation functions and the 
masses of the lowest-lying glueballs in the SU(2) gauge theory are calculated 
from uncooled and cooled configurations. The simulations are performed,.an 
an 8 3 x 24 anisotropic lattice by using tadpole-improved Symanzik's actiorEHMi. 
The lattice spacing in the temporal direction is smaller than that in the spatial 
direction with the anisotropy ratio £ = a s /a t = 3, so that the physically useful 
signals can be obtained before they would be undermined by the statistical 
fluctuations. The form of the lattice gauge action we use is expressed as 
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where f3 



g is the gauge coupling. is the plaquette operator and Rfj 



2x1 rectangular operator, u s and Ut are tadpole improvement parameters. 
u s is defined by u s —< (l/2)TrW sp > 1 ^ 4 , while ut is set to be 1 because of 
Ut = 1 — 0(£~ 2 ) (£ is much bigger than one). 
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In the continuum theory, the gluonic operators with quantum numbers + , 
0~,2 + arc defined by the field strength squared, the topological charge density, 
and the energy density, respectively, 



Os — (gG^) 2 , p — -e^pag 2 G f _ l „G pa , (2) 
OT = \(gG;„) 2 -g 2 G a 0a G a 0a . 
The correlation functions C(x) for the Euclidean separation x are defined 

as 

C r (x) = (0|Pr(x)pr(0)|0). (3) 

The glueball spectra are obtained by the spatial integral of two points 
functions of operators with definite J PC quantum numbers, such as 

C r (t) = f rff(0|Or(x)O r (0)|0) (4) 
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The lattice simulation of masses of lowest-lying hadrons is well established. 
The lattice version of Or are defined as Ell 

O 0+ (t) = E Pl2 ^ *) + P23 ( f < t) + P3l(x,t)}, (5) 

X 

2+ (t) = X! t P 12 ( f > *) - P 13 (*.*)]■ (6) 

a? 

The corresponding glueball masses are extracted from the exponential fallof of 
the lattice correlation functions. 

mr(t) = > (t^t) (?) 
a t \C(t-a t )J 

We first calculate the lowest-lying glueball masses on the uncooled configura- 
tions at (3 = 1.0, 1.1, 1.2 through the standard procedure, where the fuzzy and 
smearing methodafl are used to increase the signal-to-noise ratio. The con- 
figurations are thermalized through 5000 Monte Carlo sweeps ( each sweep is 
composed of four heatbath sweeps and one micro-canonical iteration) and the 



measurements are performed on 10000 configurations. We analyze the statis- 
tical error by dividing the 10000 configuration into 20 bins. The physical scale 
is set by calculating the heavy quark potential. With the value of the string 
tension y/a = MOMeV, a s takes the value 0.277 fm, 0.233/m, and 0.183/ttEI 
respectively for j3 = 1.0, 1.1, 1.2. The glueball masses at (3 = 1.0 are extracted 
to be 

m 0+ = 1634 ± 60MeV 

m 2 + = 2305 ± 72MeV 

m - = 2574 ±224MeV 

Our results are in good agreement with that of other worlEl. 

To explore the role of instantons in the gluonic correlation functions and 
glueball masses, the cooling method is used to get the instanton configurations. 
The link-updating of the cooling iteration is defined as below 

where c is the normalization factor so that U'^ is an element of SU(2) group. 
E M is defined by the local action 

s(J7 M ) = 1 - ^Trll^^istaples) = 1 - ^Trl^X^, 

This transformation justifies that the local action will be minimized. We take 
the measurement on a sample of 500 configurations generated and cooled by 
the lattice action S\. The gluonic correlators are measured every five cool- 
ing sweeps and the spectra are extracted by the same approach as that was 
used in the uncooled case. In order to avoid the auto-correlation, every two 
adjacent configuration are separated by 20 Monte Carlo sweeps. Each con- 
figuration is cooled by 100 steps. The cooling procedure is monitored by cal- 
culating topological charges at each cooling step. The lattice charge density 
operator(Q cont .(x) = TrF^ (x)F^ (x) in continuum theory) used in this 
work is the improved version, 

51-b (l Qp{x) l l 



Ql( x ) = -T^\oQp( x )-kQs(x)) (8) 
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Qp{x)= J2 Zw T r{PA x ) p P*(. x )) ( 9 ) 
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Figure 1: The correlation function C(t) in 0+ channel at various cooling steps arc illustrated. 
The dots, the rectangulars, the triangles and the upside-down triangles correspond to the 
cooling steps 10, 25, 50 and 100, respectively. 



Qs(x) = £ WTr [(S^ 2 (x) + Sl*\x)) (S% 2 (x) + S^(x))} . 
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Here e^po- is the standard Levi-Civita tensor for positive directions while for 
negative ones the relation e^ vpa = —t-^vpa holds. is the plaquette opera- 
tor in the /j, — v plane, while Sffi is i x j rectangular operator. We find that 
after several cooling sweeps the charges will reach meta-stable plateaus with 
approximate integer values. 

Fig. 1, 2, 3 show the correlation functions of the three channels. In each 
channel, the correlation function is plotted at 10, 25, 50 and 100 cooling steps. 
The glueball "masses" are extracted in the same way as in the uncooled case 
and are plotted in Fig. 4, 5, 6. 

Shuryak and Schafer has calculated the gluonic correlation functions in 
the single-instanton background field. TakingJthe free gluon propagator into 
consideration, the correlators are expressed asU, 

C s (x) = ^4 + np*C inat {x), (10) 
C P (x) = + np 4 C mst (x), (11) 

Ct(x) = (12) 



C 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
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Figure 2: The correlation functions C(t) in 2+ channel at various cooling steps are illus- 
trated. The dots, the rectangulars, the triangles and the upside-down triangles correspond 
to the cooling steps 10, 25, 50 and 100, respectively. 

where g is the running coupling constant and Ci ns t {x) is the instanton contri- 
bution. Ci ns t (x) is definitely positive and its detailed expression is abbreviated 
here. The signs in Eqn.(10) and Eqn.(ll) are impressive: the two contribu- 
tions reinforce in the scalar channel(attractive) but cancel to some extent in 
the pseudoscalar channel(repulsive). There is no effect in the tensor channel. 

Our results shown in the figures are qualitatively in agreement with these 
comments. At small t, the amplitudes of the correlation functions decrease 
rapidly with the increase of the cooling steps, so do the extracted masses. In 
the intermediate and large t range, the correlation functions are less sensitive 
to the cooling. The "masses" take out short plateaus in the + channel and 
the 2 + channel in the intermediate t. These are the results of that the ul- 
traviolet components of the gauge fields are removed by the cooling but their 
topological structure is kept more stable. Small instantons are smoothed out 
by the cooling iterations (we monitor the cooling procedure and find that only 
the instantons with sizes greater than 2a s survive), so the correlation functions 
and the "masses" at small t may be underestimated. The fact that in the 2 + 
channel the correlation function tends to be flat (Fig. 2) and the mass-platcau 
decreases rapidly (Fig. 5) imply that there is no instanton effect in this chan- 
nel. In contrast to this, the fall-off the correlation functions with t in the + 
and CP channels are obvious even after 100 cooling sweeps. The height of the 
mass plateau in the + channel is less sensitive to the cooling steps in the 
intermediate t range(Fig. 4). This means the contribution of instantons to 
the glueball mass in the + channel is significant. As illustrated in Fig. 3, the 
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Figure 3: The correlation functions C(t) in - channel at various cooling steps are illus- 
trated. The dots, the rectangulars, the triangles and the upside-down triangles correspond 
to the cooling steps 10, 25, 50 and 100, respectively. 



correlation function in the 0~ channel even changes sign at t = 5at — Sat after 
50 cooling steps. This is an evidence of that the gluonic interaction induced 
by instantons is repulsive in the pseudoscalar channel. The mass information 
can not be clearly extracted in this channel (Fig. 6). 

3 Conclusion 

We simulated the SU(2) glueball masses in the cooled and cooled gauge config- 
urations. The results in the uncooled case are in agreement with the existing 
work. In the cooled case there is evidence of the instanton-induced attractive 
and repulsive interactions in the + channel and the channel, respectively, 
which are predicted by Shuryak and Schafer based on the instanton liquid 
models. The contribution of instantons to the + glueball mass is significant. 
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Figure 4: The 0+ channel 'glueball masses' are extracted from the exponential decay of 
the correlators in cooled case. The masses vs. Euclidean time at various cooling steps are 
plotted. 
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Figure 5: The 2+ channel 'glueball masses' are extracted from the exponential decay of 
the correlators in cooled case. The masses vs. Euclidean time at various cooling steps are 
plotted. 
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Figure 6: The 0~ channel 'glueball masses' are extracted from the exponential decay of 
the correlators in cooled case. The masses vs. Euclidean time at various cooling steps are 
plotted. No plateaus appear. The mass information can not be obtained in this channel. 
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